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Work, under NASA Grant #NAG 1-158 has progressed signifi- 
cantly during the past year. The original purpose of this 
grant program was to investigate the nonlinear and dynamic 
response of composite structures to impact loading conditions. 
To understand this problem, it was necessary to investigate the 
fundamental properties of the deformation processes. To this 
end, the work initially focused on the response of isotropic 
beams subjected to pulse loading. The convergence of both 
finite element and finite difference solution procedures was 
investigated with this configuration. 

Existing solutions in the literature did not model the 
experimentally observed behavior of impacted plates. The 
convergence studies demonstrated that the reason was due to the 
lack of an appropriate error norm. The total system energy was 
investigated and determined to be a good measure of accuracy. 
This norm was used for all further studies. The work on 
isotropic beam systems resulted in two publications listed as 
Appendixes A and B of this report. 

Experimental studies were also conducted to examine the 
response of composite beams and plates to impact loading. The 
experimental work was conducted during the summer of 1984 at 
the Langley Research Facility under the joint direction of 
Dr. Wolf Elber (of NASA Langley) and Dr. E. Thomas Moyer, Jr. 
(of The George Washington University). These studies were 
carried out by Ms. Ingeresa Pincus. Here tests provided 
fundamental data on the indentation and flexural response of 
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the materials being studied in both beam and plate geome- 
tries. In addition, several fundamental properties of the 
deformation and onset of failure were elicited by these 
studies. The results of her work, are presented in Appendi* C. 

The numerical analyses were extended to the study of 
circular plates subjected to impact loading. The preliminary 
studies involved a critical computational investigation into 
the behavior of the finite element method and finite difference 
method for solving dynamic nonlinear flexural response 
problems. The results indicate that the finite difference 
method (though rarely used in solid mechanics anymore) can be 
very efficient for solving nonsingular problems provided a 
vector computer is available (e.g., the CYBER 205). Finite 
element solutions also were generated. For the problems 
studied, both methods were comparably efficient, however, due 
to the simplicity of the finite difference method and the easy 
vectorization, it was preferable. The comparison studies were 
the topic of an invited conference paper included as Appendix D 

Further work was done in modeling the response to simu- 
lated impact loading through careful modeling of the momentum 
transfer from impactor to target. The extreme computation 
times required for these solutions limited the example problems 
which could be solved. The modeling methodology and solution 
procedures were, however, well tested and verified. The re- 
sults of this work resulted in an invited paper presented at 
the ASME Winter Annual Meeting (1984) and is included as 
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Appendix E of this report. 

To attempt to reduce the computational requirements, 
fundamental investigations into the convergence, stability and 
accuracy of several common time integrators for dynamical 
systems was initiated. The results demonstrated that the 
Runge-Kutta integrators of order four may prove to be more 
competitive for flexure problems than those commonly employed 
(e.g., central difference method, Newmark's method, etc.). 
Preliminary results were presented at the SES Annual Meeting 
(October, 1985) and will be the topic of a future NASA 
contractor report and publication. The results presented at 
the meeting are summarized in Appendix F. 

During the summer of 1985, the focus of the research 
changed to investigate the response of plates to acoustic 
loading. Preliminary unpublished experiments from the NASA 
Langley Reseach Center demonstrated that the response of many 
panels relevant to application exhibited significant nonlinear 
characteristics. The solution methodology developed for the 
impact simulations provides the basic tools necessary to 
address these problems. Our most recent work, is focusing on 
the development of accurate and efficient solutions to acoustic 
response of plates. The formulation of the fundamental equa- 
tions and the specific problems being investigated are pre- 
sented fully in the project proposal. 

/ 

A computer program for the solution of circular plate 
problems has been developed and is under testing. The program 
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can analyze arbitrary loading and damping properties. Finite 
difference operators have been employed for spatial discreti- 
zation. Different time integrators are being utilized as 
necessary as the damping terms introduce different convergence 
properties than exhibited by the undamped system. 

In addition to the work, cited above, further research has 
been done in the convergence properties of finite elements for 
solving many types of nonlinear vibration problems. This work 
has delineated the grid requirements, time step requirements 
and error norm magnitudes necessary for accurate solutions to 
be generated. This work was performed by a candidate 
(Mr. E. Ghashghai-Abdi ) for the Master of Science degree at The 
George Washington University and is the subject of his thesis 
(in progress). This work will be the subject of a NASA 
contractor report at its completion. 

The work under NASA Grant NAG #1-158 has produced many 
important results summarized above and detailed in the 
Appendices. In addition, the computational approaches devel- 
oped for solving impact problems are largely applicable to 
acoustic problems. The extension into this domain was natural 
and is the subject of our current work. The investigators wish 
to thank Dr. Wolf Elber (Technical Monitor of this Grant) for 
his continual support, encouragement and involvement in this 


research. 
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SUMMARY 


Static indentation tests were performed on 8-ply quasi-isotropic 
graphite /epoxy composite materials to serve as a basis of comparison to dy- 
namic tests. 

Low velocity impact tests were conducted with the same materials 
to provide supporting data for theoretical analysis developed by Dr. Thomas 
Moyer. Strain was recorded for several locations on both circular and 
beam-shaped specimens. Strain versus time graphs were then plotted for 
impacts of various magnitude. 

The results demonstrated that the strain behavior was a function 
of position along the radius and dependent upon the flexure and membrane 
behavior of the material. 


INTRODUCTION 

The increasing use of composite materials in aircraft structures ne- 
cessitates research into deformation and failure patterns. Determining dam- 
age from low velocity impacts is of great importance because the occurrence 
of this type of impact is highly possible. In addition, resulting damage 
which could propagate through the material may not be visible on the sur- 
face. 

The purpose of this report is to explain the static indentation and 
low-velocity impact tests conducted on composite materials. The report will 
include a description of the test set-ups and procedures and and a graphical 
presentation of the test results. 



EXPERIMENTAL PROCEDURES 


Materials 

Two 8-ply quasi-isotropic composite materials of graphite in Narmco 
5208 epoxy resin (T300/5208) were tested. Material A was a woven fabric 
laminate and Material B was a tape laminate. Specific material properties are 
listed in Table 1. 

Static Tests 

Static indentation tests were conducted first. These tests, de- 
scribed briefly below, generated load versus displacement curves and strain 
versus displacement curves. 

Test Set-up 

The composite materials were cut into 7" squares and l"x7" beams 
and bonded with epoxy to square aluminum plates with 5 n diameter aper- 
tures. The specimens were clamped to a test stand mounted onto the hy- 
draulic ram of an MTS servo-hydraulic test machine. The indenter, a cylin- 
der with a 1" steel ball attached to the tip was mounted to a load cell from 
above. Two DCDT displacement gauges attached to the sides of the indenter 
measured the displacement of the plate under loading. This set-up is shown 
in Photos Nos. T and 2. 

Test Procedure and Analysis 

During one static indentation test cycle, the displacement followed 
a triangular waveform with a set maximum. The resulting load applied to the 
specimen during the test was recorded. (The calibration of the displacement 
gauges was 1 volt = 0.6312 mm. and the span of the machine ranged from 0 
to 10 volts.) The load versus displacement curves generated for the two 
materials are shown in Figures 1 and 2. 
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According to the analysis developed by Dr. Wolf Elber and 
Gretchen Bostaph in Ref. 1, the reaction of a plate to a load is determined 
by the combination of plate shear and flexural stiffness, together with mem- 
brane extensional stiffness. The theory predicts that large diameter plates, 
such as the test specimens used, would be dominated by the membrane be- 
havior. The analytical load versus displacement curve for a mem- 
brane-dominated plate, presented in Ref. 1, has an initial continuous, non- 
linear slope as did the curves in Figs. 1 and 2. Comparatively, the typical 
curve for a small diameter plate, where the flexural behavior dominates, has 
an initial linear slope followed by a sudden drop, indicating the onset of de- 
lamination . 

Strain gauges were mounted at several locations along the radius 

of the plate and measured during small displacements of the plate. Strains 

were plotted against displacement, demonstrating changes due to flexural and 

membrane strains. Figures 3 to 8 demonstrate strains on a circular plate 

and Figures 7 and 8 show the strains on a fixed beam. As predicted in 

Ref. 1, the flexural strains were tensile inside the inflection point and 

compressive outside, for the bottom surface. Strains on the top surface were 

/ 

opposite in sign. The overall strain behavior, however, was also influenced 
by membrane tension. The dominating contribution of membrane strains was 
exhibited by the graphite epoxy material, shown in Figure 5. 

The indentation of the plate, caused by the indenter, contributes 
to the total displacement of the plate. As explained in Ref. 1, a function 
representing this indentation is determined from the load versus displacement 
curves of specimens totally supported underneath. Unlike the other static 
indentation tests, the machine was set on load control and the displacement 
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recorded. (One volt corresponded to 400 pounds of load.) These indenta- 
tion functions (representing the load versus displacement curves shown in 
Figures 9 and 10) were calculated by transforming the data into linearized 
form and performing linear regression. The linearized data and the best line 
fit is shown in the Appendix. 

Low Velocity Impact Tests 


Test Specimens 

Two types of specimens were tested: 5 n diameter circular plates, 

similar to the status test specimens, and beam plates l"x5" , fixed on the 
ends. Eight strain gauges were bonded to each specimen in specific lo- 
cations with respect to the point of inflection. Strain is initially zero at the 
inflection point for a plate subject to a center unit load. The exact location 
was calculated from the following equation for deflection of a circular plate 
with clamped edges, loaded at the center: 

i.r _ Pr 2 log r/a ^ P(a 2 -r 2 ) 

W 8ttD + B'uD ' 

where "a" is the radius of the plate, and "r" is the distance from the cen- 
ter . 

The second derivative with respect to "r" is: 

d^w _ P [1+log r/a] 
dr 2 4 n D 

From, this equation the location of the inflection point was ob- 
tained. 

lnr = -1 + lna 
r/a = e * = 0. 37 

As shown, the point of inflection was located at a distance 
r = 0.37a from the center. Strain readings were desired at this point in or- 
der to record the fluctuations in strains due to the combined effect of the 
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membrane and flexural strains. Strains were also recorded at a distance 
r = 0.75a where the strain behavior was less sensitive to position. In order 
to observe differences between strains oriented along the direction of the 
fibers (0° direction) and strains oriented at 90° from the direction of the 
fibers, gauges were placed in both directions on the circular plates. The 
exact locations and labels for the strain gauges bonded on the 127 mm. cir- 
cular plate were the following: 

- Strains A & B: r = 23.5 mm., top and bottom surface 

- Strains C & D: r = 47.6 mm., top and bottom surface 

Strains A through D were all oriented in the 0° direction on the 
circular plates and in the 45° direction on the beam plates. 

Strains 1 through 4 corresponded to the same locations as A 
through D, except that they were oriented in the 90° direction. 

Test Procedure 

The device used to produce the impact consisted of a l" diameter 
steel ball attached to the end of a fiberglass cantilever. Mounted on top of 
the ball was an accelerometer to measure the vibration of the plate. The 
specimen was placed in a position so that the cantilever rested horizontally 
and unstressed when the ball contacted the center of the specimen. This 
test set-up can be seen in Photos Nos. 3 and 4. To produce the impact, 
the cantilever and ball were lifted to a specific height and dropped. The 
release heights covered the range of 1" to 8" in one-inch increments. Each 
height corresponded to a specific maximum load to which the plate was sub- 
jected. The equation for this "max” load was derived by representing the 
cantilever as a simple harmonic spring and determining the acceleration in 
terms of the duration of the impact. This max load is: 


Pmax = 


-T> VoM 

T 
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Where "M n is the mass of the ball and cantilever, n T n is the duration of the 
impact, and n Vo" is the velocity of the ball when it strikes the plate, which 
is represented by the equation: 

Vo -VP + 2 * h 

where "K" is the stiffness constant for the cantilever and n h" is the release 
height. The constants used for analysis are found in Table 1. 

DATA 

Strains were read and recorded on a four channel Nicolet Oscillo- 
scope. The strain gauges used had a gauge factor of 2.010.1 and a resis- 
tance of 120 ohms. The output was measured and amplified through a 
Bridge Amplifier and Meter (BAM-1) and calibrated for 1 volt equal to 
4000/? £. These voltage readings were recorded on Channels 1A and 2A of 
the oscilloscope. The time duration of the impact was recorded on Channel 
2B and the acceleration on Channel IB. The equipment used can be seen on 
Photo No. 5. 


ANALYSIS 


The information recorded on the four channels of the oscilloscope 
was originally stored on a magnetic disk. During the test, readings were 
taken every two milliseconds, producing an extremely large number of data 
points. The data was fed through a program designed to reduce the number 
of data points to below 200. The new data set was then stored on magnetic 
tapes. The strain readings were then converted from volts to microstrains 
by a calibration program. Strain versus time curves were plotted to display 
the test results in useful form. 


- 6 - 



RESULTS 


The results of the low velocity impact tests are shown in Figures 
11 through 16. The strains measured for each impact were plotted according 
to the location of the strain measurement. The different curves represent 
the increasing release heights for the impacts, producing increasing loads on 
the specimens. 

The graphs indicated that the duration of impact decreased as the 
"max" load of impact increased. The larger impacts caused an increase in 
strain magnitude. However, the pattern of strain behavior remained very 
similar in most cases. The large contribution of membrane tension in the 
graphite epoxy composite is illustrated in the dynamic test curves by the 
tendency towards positive total strain at locations of compressive flexural 
strain . 

The results of the beam specimen differ from those of the circular 
plate specimens primarily in the magnitude of the strains. The beams re- 
spond with higher strains and also illustrate more flexural characteristics. 

CONCLUDING REMARKS 


An 8-ply graphite epoxy composite material and an 8-ply woven 
fabric composite material were chosen for testing. The materials were cut 
into 5" circular plates and S'^l" beams. 

The following tests were conducted: 

1. Static indentation: From these tests, load versus displace- 

ment curves and strain versus displacement curves were gen- 
erated. 
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Indentation on fully backed specimens: These tests produced 

the indentation functions for the two materials. 

Low-velocity impact: These tests provided the strain behav- 

ior data needed to support theoretical analysis developed by 
Dr. Thomas Moyer. 
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